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Abstract 

Here, we extend the analysis of part I to determine the linear stability of a bubble in 
a Hele-Shaw cell analytically. Only the solution branch corresponding to largest possible 
bubble velocity U for given surface tension is found to be stable, while all the others are 
unstable, in accordance with earlier numerical results. 
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1. Mathematical formulation and regular perturbation for the modes 

Now we consider the problem of linear stability of a bubble for disturbances that 
are symmetric about the channel centerline. For the class of disturbances considered, 
the flux of fluid at infinity is assumed to be the same as in the steady state. This is a 
realistic assumption for most experiments where the pressure gradient is held fixed. From 
incompressibility of fluid flow, the bubble area is time independent. Our formulation here 
is an extension of earlier method for the analytical determination of the zero surface tension 
symmetric stability modes as shown in section III of Tanveer & Saffman(1987). 

Without any loss of generality we set the channel half width a to unity and the fluid 
velocity at oo , V to unity since this is equivalent to nondimensionalization of all relevant 
variables using combinations of these two quantities. In the frame of the steady bubble 
translating with constant velocity U , the pressure boundary condition on the bubble 
boundary is 


( 1 ) 


b 2 T 1 

6 + Ux = — + constant 

12 fiR 


where the variables have the same physical meaning as the first part of the paper. The 
constant appearing on the right hand side of (l) can now depend on time. The kinematic 
boundary condition that a point on the bubble boundary moves along with the velocity of 
the fluid at that point can be expressed as 


(2) <f> x F x + <f> v Fy + F t = 0 on F(x, y,t) = 0 

where -F(x,y,t) = 0 determines the bubble boundary (x(t),y(t)) at time t . As for 
the steady problem (Tanveer, 1986), we consider the mapping function z($, t) from the 
interior of the unit f semi-circle into the fluid flow region exterior of the bubble on one 
side of the channel centerline (z = x + iy ) such that f = ±1 are mapped to the rear and 
front stagnation points on the bubble respectively; f = ±a are mapped to z = ^oo . 
The semi-circular arc corresponds to the bubble boundary, while the real diameter in the 
f -plane corresponds to the wall y = 1 and the channel centerline y = 0 . a depends 
on time and has to be determined. In the linear stability problem, it is convenient to 
decompose 

(3) a = a + <*! 

where a is time independent and is specified as a parameter characterizing the size of the 
steady bubble whose stability is being studied. It is convenient to decompose the mapping 
function into 

(4) 2 = z 0 ($, a) + t) 
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where f is some analytic function within the unit circle that has to be determined and 
zq is given by 

(5) *o(f, «) = (^f) + J (| - l) ** (£ff) 

and k t is a constant given by 

( 6 ) 


2a 


We find 

( 7 ) 

where 

( 8 ) 


k. = ^lU(l + a>)-2a>] 


= *.(/' + r)h ) 


= a[U{ 1 + a 2 ) ~ 2a 2 ] 
** a[U{ 1 + a 2 ) - 2a 2 ] 


and the function h, is given by 
(9) h = ^ ^ 


where p 2 = ^ + *> ' “ 2 


(f 2 - a 2 )(l - a 2 ? 2 ) ' 17(1 + a 2 ) -2a 2 

We decompose the complex velocity potential W = <f> + i 0 into 

(10) W = W 0 (f,a) + C//: 4 u;(f,t) 
where 

( 11 ) 


Wo(f,a) = — UfalL “M 1 £!£ 

v 7 T (f + a)(l + af) 


On substituting (10) and (5) into the dynamic boundary condition (1) on the bubble, we 
find that it is equivalent to 

— / 1 -t- Rzf — 

( 12 ) 


„ , , {1 + Resf tn(!‘ + rih)} 

He(W + /) = - 1 


on f = e iv where u is real and in the interval [0, 7r] and 

b 2 T 

< 13 > 

as in the first part of the paper. Since for symmetric disturbances, the streamlines coincide 
both with the cell walls and the channel centerline, 


(14) 


Im u> = 0 
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on the real f diameter. The geometric condition that real f axis in the interval 
[—1,1] corresponds to the wall and the channel centerline is equivalnet to 

(15) Im f = 0 


on the real f axis in [-1, 1]. Since |$j = 1 corresponds to the bubble boundary at all 
times, the kinematic boundary condition (2) is equivalent to the following equation on the 
unit £ semi-circular arc (see Tanveer & Saffman (1987) for details) 



where * indicates complex conjugate and subscript t denotes time derivative. On 
substituting for z from (4) and W from (10) and using (ll), one finds that (16) is 
equivalent to the following boundary condition on f = e w u in the interval [0, tt] : 

(17) Re[ f w'-f(/' + oi)-(g + /),]= 0 

where in the above prime on quantitites u and / denotes derivatives with respect to 
f and subscript r indicates derivative with respect to r , where r is scaled time given 

by 


(18) 


T — 


Ut 

k. 


For linear stability analysis, we decompose 


(19) 


/ = /+*• 


where / satisfies the steady state equations (1) and (2) in the first part of the paper. 
Assuming F and a\ small compared to / and a , linearization of (12) and (17) about 
the steady state gives us the following boundary conditions on f = e w : 

( 20 ) 


Re (w-hF) = 


\f' + h\ {Re *d< 


d + _ 


r+h j 


Re 


1 + ^ In ( r + h ) 


Re 


F' + Hiai 


f' + h 


} 


( 21 ) 


Re 


+ F,) - i, + F t ) 


‘ ' ' k. 

where subscript a denotes derivative with respect to a . 


0 


( 22 ) 


Mf) = 


i -pY 


(f 2 - a 2 )(l - a 2 f 2 ) 


4 


_ (1 — p 2 i 2 ) 1 2aU(l — f 2 ) — 4a 2 a 2 a$ 2 

1 (f 2 — a 2 )(l — a 2 f 2 ) a U( 1 + Of 2 )(l — f 2 ) — 2(a 2 — f 2 ) £ 2 — a 2 **" 1 — a 2 f 2 

For f real in (-1, l), in addition to (14) substitution of (19) into (15) and use of equation 
(l) in Part I implies 

(24) ImF = 0 

Equations (21), (22), (14) and (24) are the linear stability equations determining the 
analytic functions u; and F and the time dependent constant ai . For purposes 
of analysis of these equations, it is convenient to define analytic functions k , R and 
T within the unit circle with continuity at the boundary such that on the real diameter 

(25) Im k = 0 


Im R = 0 


Im T = 0 


and on £ 
(28) 

Re k = - 


r, +'>i{ 


d ( F f + H\cci \ d . . N F 

Re [ TTh ) ~ Re [ 1 + ln (/ + h) Re ~ 


f' + h JJ 


ReR = Re k + 


(30) Ref = Re[?f'*F] 

Each of k , R , T so defined are analytic functions of £ in |f| < 1 because of smoothness 
of the right hand side of (28), (29) and (30) on f = e xv and from the application of 
Schwarz reflection principle for continuation to the lower half f semi-circle. Substituting 
for h from (22) into (20) and using the definitions of k , R and T it follows that on 


<r = e* 


Re [oj + F — qr/c] = 0 


Re kw' 


(Jl ~ P 2 ) 

(< T 2 “ « 2 )(1 ~ « k 2 ) 
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By subtracting off the simples poles of the expression within the square parantheses in 
(32) and adding on regular terms with the same real part, one finds that (32) is equivalent 
to 


(33) 


Re S = 0 


where 

(34) 


s = 'u> | F + p„, +T | i 

C + (l-«y)({*-a*) F ' +fi ‘' +T,+ 2 (l-a«) (f 


*V(«) 


F,(-a) 


(? - a)(l - af) (f + a)(l + af)J 


Clearly Im S = 0 for f real in [-1, 1] because of (14), (24), (26) and (27). Therefore, 
(33) implies that S = 0 for |f| < 1 and so 
(35) 


(k 2 ~P 2 ) r, , f,. , * . ~ P‘) ,.t 

f + (l-a*(’)((*-a ! ) F,+R u+Tr+ 2(l-a<)k 


*V(“) 


F r (-a) 


(C - a)(l - as) (s + a)( 1 + as) 

Again the function within the parantheses on the left hand side of (31) is an anlaytic 
fucntion in |f| < 0 with vanishing imaginary part on real f axis in [-1,1] in view of 
(14), (24) and (25). Thus for any f with |f| < 1 , 


= 0 


(36) 


u + F — 7 /c = 0 


Equations (35) and (36) hold for any f with |f| > 1 as well when one obtains proper 
analytic continuation of each of k , R and T defined originally in |f| < 1 . On 
substituting for u> from (36) into (35) one obtains 
(37) 

■ r.. . S fc ~ P 2 ) p .p ,f . ( a, -P‘) l .2 *V( a ) . F,(-a) 1 _ - 

K * - a)(l - as) + (f + «)(1 + «r) j ~ 1SIC 


On evaluation of the above equation at f = 0 , we find that 
(38) 

and this determines O'! . One can use (38) back into (37) to obtain 

rr-^V = Rs r + IK 1 


2 ( “ ( , _^] ) [-FV(a) - F,(-o)] + T,( 0) + X(0) air = 0 


(39) 

F' + ^ 

(1 - a* 

where 


(40) 

p n<*)h 

Jt3 — 


(i ! ~ P 2 ) 


+ 


? + a 


+ 


i2(0)f 
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where 


(41) 


(42) 


P _ (* 2 ~ P 2 ) 
1 2(1 - a 4 ) 


C 2 ~ 1 + (f - a)-R(f) 


A = 


(<* 2 - p 2 ) 
2(1 - a 4 ) 


[ <r(l -a<r) ' a R(0)c 

f 2 - 1 (f + oc)R{s) 


f(l + <*?) a£(0)f 

We note that .R(O) ^ 0 is a necessary condition for the above to be valid. This conditions 
is satisfied when 7 = 0 for in that case R = R 0 , where from (5), (22) and (29), Ro is 
determined from the boundary condition 




f 2 ) 


+ 


C(« 2 - P 2 ) 


[ f 2 — a 2 (1 - a 2 p 2 )(l - a 2 f 2 )J 


1 


on f = t xv and on the real diameter the imaginary part of R 0 vanishes. It is not 
difficult to see that 

1 («* - P f) 


Ms) = 4( ? , 

7r k 8 { (f 2 - a 2 )(l - a 
, (<* 2 - p 2 ) 


- a 2 f 2 ) 


+ 


( C 2_ a 2) (l- a 2^2)( 1 _ a 2 p 2)J 


(44) 


4(1 - a 4 ) 


+ 


L(? -a ) 2 (f + °0 2 (l-<*?) 2 (1 + af) 2 


^ / 1 2a ^ 2 a s \ ( 1 f 1 $■ \ "I 

\a a 2 — p 2 1 — a 4 /\f — a 1 — af f + a 1 + af/J 

i_- -U) 

fa 1 + afJJ 


a(a 2 - p 2 ) 


-L.+-J-+ 

C-a 1-af f + 


and so 

(45) 


2(1 — a 4 ) 2 (l — a 2 pi 2 ) 

R( 0 ) = — [ ( a2 ~ p2 ) ( 1 + . (a 2 ~ P 2 ) 2 ) 

0 ttA:, \ (1 — a 4 ) \a 2 — p 2 1 — a 4 y (1 — a 4 ) 2 (l — a 2 p 2 ) J 


which is clearly non zero. For small 7 , / is small and therefore one could argue that 
i2(0) ^ 0 for sufficiciently small 7 . We will assume that is indeed the case. Now we 
assume that F has a e at time dependence then the mode corresponding to a , we can 
try a regular perturbation expansion of the form 


(46) F = F 0 ($ } e,U,a) + 7Fi(<r,a,l7,a), + ... 

where for convenience of regular perturbation expansion, U is being treated as a parameter 
together with a , though the steady state selection implies that U is a function of 7 and 
a . Tanveer & Saffman (1987) calculated the eigenmode F 0 . Using their procedure, it is 
straighforward though laborious to find the next term F t and one finds no restriction on 
a . We assume that this is indeed the case to every order in 7 . 
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2. Analytical continuation and reduction to an ODE 

We now proceed to extract the leading order transcendent ally small terms by analytic 
continuation of the boundary conditions to the unphysical region |f | > 1 . The procedure 
is in principle very much like the first part of this paper. This illustrates the power of this 
formalism, which otherwise may appear rather complicated compared to the Combescot- 
Dombre analysis. The complication here compared to the first part of this paper is that 
we have two analytic functions / and w rather than one. We now proceed with the 
analytic continuation of the equations. Rather than work on the linearized set of equations 
just derived, it is convenient to analytically continue the full nonlinear equation and carry 
out the linearization at a later stage. 

We can write equation (12) as 


(47) 


Re 


! | s :(r + jfr) . g 

(/' + rih) 7. 


= 0 


where g is an analytic function in the unit semicircle defined by the boundary conditions 


(48) 


Re g — Re(w + f) \f' 4- r)h\ on ( = p e w 


and on the real diameter, 


(49) 


Im <7 = 0 


Using Poisson’s integral formula relating a harmonic function and its conjugate to the 
value of the harmonic function on it’s boundary, by using the definition of h , we find that 
(48) implies 
(50) 





(f ' 2 - a 2 )(l - &V 2 ) 



+ w(f') + u 



where 


(51) 2i(f) = xfV - p j ) + (f ! - S ! )(1 - &V)/'(i) 

and 


(52) U() = 1(1 - lV) + (f 1 - &’)(! - &Y )'m 


Analytic continuation outside the unit circle gives us 


(53) 


g{c) - |/(f) + f(-) + w (c) + w(-)| 
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Now, we notice that the quantity within the paranthesis in (47) is an analytic function of 
f with simple poles at £ = ±a and has vanishing imaginary parts on the real diameter. 
By subtracting off the simple poles and adding on terms with the same real parts on the 
unit circle, it is easy to conclude that (47) is equivalent to 


(54) 


f(7* + lfr) . 9 . 

(/' + r)h) 1 f 2 - d 2 


2 

1 - d 2 f 2 


= 0 


Equation (54) is valid everywhere in the complex plane. For | > 1 , we must use 

the expression (53) rather than (50) for g . Using the decomposition (3) and (19) and 
assuming that each of ai , F and u> are small, we expand I to linear order to find 


(55) 
where 

(56) 

where 

(57) 


/(/,«,£) = Hf,() + + «,/„(/, f) 


1 1 


f - cj b" - £>’)(! - a*(”) (/(f 1 + 


<i(f) = f 2 (f 2 - P 2 ) + (f 2 - ° 2 )(1 - aV)/'(i) 


(58) «,(?) = (1 - pY) + (f 2 - a 2 )(l - a 2 fV(?) 

and 


(59) 


{, hl f 


<r + <r' 


?'-?j 


{*•(/.**) 


f( f 0 + w ( f ') + f(p) + «(i) 


+ ir < (/,? , )F , (f’) + /f 5 (/,f')i 2 '(p)} 


where 


(60) 

(<r 2 -a 2 )(l-aV) 

(61) 

W,?) = ^k i [/(C) + /(i)l 

(62) 

«i(/,f) = |<r*<l[/(f) /(|)i 
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and 


(63) 


r „ , i /d<r'ff + f'l4k')4(f , )(/(f') + /(?)) 




(?' 3 - a 2 )(l - a 2 ? 33 ) 




2a 2a?' £i(?',a) £i(?',a) 

4 . »2 I / .i\ ~T 


where 


^_ a 2 i_ a 2^ 2£ 1 (<r') 2£ 2 (f') 




(64) 


" fl { 


2Ua 


2a(U(l + a 2 ) — 2)(£7 — 2) 


£7(1 + a 2 ) - 2a 2 (17(1 + a 2 ) - 2a 2 ) 2 


-2a(l - aV)/'(i) - 2a?V - a 2 )/'(i) 


£ i(? 




2a(£7 - 2) 
£7(1 + a 2 ) — 2a 2 


}d-p J 


f ! ) 


(65) 


+ 

2 / 2Ua 2a(U(l + a 2 ) - 2) (£7 - 2) ) 

C \ £7(1 + a 2 ) -2a 2 (£7(1 + a 2 ) - 2a 2 ) 2 J 


-2a(l - aV)/'(f) - 2 af 2 (f 2 - a 2 )/'(f) 

Linearizing (54) and using (53) and subtracting off the steady state, we find 

( 66 ) 

where 


F" + (P, + Q,)F' + ^(F + u>) = Rz 


(67) Ri = a.P.fe) - <Z!±^1/ 5 (P, w ,/,?) - ^{P(i) + „(!)} - I(71±^ Xs (/, { )p'(I) 


where 

(68) 


-(/” + £') 1 2? 
-Tl = — . , V = - + 


2 0 
+ — 


(/' + ^) £ ? 2 -a 2 f(l - a 2 f 2 ) 7^ 

where g is as defined in the first part of the paper and so from equation (20) of Part I, 


(69) 


9_ UJjSI + ry^) + ft—)] 

IS IS 1 S(S 2 ~cc 2 )( l-a 2 f 2 ) l/ig + / V JJ 
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and 


(70) 
and 

(71) 
and 


< 3 . = + f( 1 -)) 


Ll = l?+V 


44 


B, — 


f (f ! - a ! )(l - «V) 


A ( » \ _ *<>(' 

s di\f' + h) (f ! - o 


4ccf 2 


: 2 ) 2 (1-aV) 2 


(72) 


mm 4(t)4(t)m + /(f)) 

1 l(f S -<* ! )(l-aV) 


f 2 a 2af 

{ f 2 — ofl + 1 — a 


Lx(C,a) , U{$ 


+ — — — + 

2 f 2 2 l x 212 


^}] 


The linearized kinematic boundary condition is given by (35). We need the analytic con- 
tinuation of each of T and R outside the unit circle. For T defined by (27) and (30), 
it follows from Posson’s integral formula that for f inside the unit circle 

(74) f =hL,y{^ 7 ) 

For | > 1 , by contour deformation in the f' plane, we obtain the analytic continuation 
of (74) to be 


(75) 


T = h{f,F,$) + 


i/'(^k) + 7/'b)f(i) 


We noted before that R = Rq given by (44) when 7 = 0, i.e. / = 0 . So from (26) 
and (29) and Poisson integral formula it follows that for |f| < 1 , 


(76) 

where 


R = h + Mf,F ti ) 


( 77 > ^ f (£rf) [p/'(l/fVa.M + WMM 

Analytical continuation for |f | > 1 gives us 

(78) R = R 0 + I 4 (f,F,t) + ^/'(l/f)*o a (<r) +f/'(?)«b a (l/?)) 
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v- - 


Since equations (35) and (66) do not involve time r explicitly, we now assume a modal 
form with a t aT time dependence. We note that in (66) the right hand side Ri involves 
values of F and w either on the boundary of the unit circle or inside where correction 
to the regular perturbation expansion 


(79) u 0 ($,a,U,a) + 7Wi($,o,U,a) + ••• 

for u> and the corresponding series for F (given by (46)) for F are transcendentally small. 
Thus one can replace F and u on the right hand side of (66) by its regular perturbation 
expansion. In our case, as far as obtaining leading order transcendental correction, it is 
only necessary to retain the zeroth order approximation F = F 0 , u = w 0 in evaluating 
J?2 • Similarly we can substitute the regular perturbation expansion for F and u in the 
integral term Is and J 4 (/, F, f) as well as for F(l/c) and R given in (75) through 
(78). On carrying out this substitutions, (35) and (66) becomes a set of two ordinary 
differential equations for unknown modes F and u; : Strictly speaking, we should use 
different notation for the modes since F and u> were originally used to denote the general 
time dependent deviation from the steady state. However, from this point onward, we will 
only be talking about the modes and so in the interest of non- proliferation of symbols, we 
retain the same symbols for the modes. £From the two differential equations thus obtained 
from (35) and (66), we eliminate the variable F to obtain a single third order equation 
for w . 

, . ... ., ^ 2 S'\ fS" S' LA . oU 

(80) u> f Pi + Qi + + ~g (-^l + Qi) + “J w v = R* 

where 

_ il ~ 

(? 2 - 

and 

/ \ — Rit nii ,2 S' . _j / S" 5'. - . L\ \ 

(82) Ri = — R Zt - (— + Pi + Ql)Rsr ~ y-g- + ^-(Pl + Ql) + — J ^3r 

where R 2 , P3 , Pi , Qi and L\ are given by (67), (40), (68), (70) and (71) respectively. 




12 


3 . Asymptotic solutions for small 7 

For small 7 , dominant balance involving terms with factor of - must give u = w 0 . 
Carrying this out this out to higher and higher precision only generates the perturbation 
expansion (79). It is clear that transcendental terms in surface tension must be generated 
by looking at the homogeneous part of (80) and trying a WKB approach. We find that 
the leading order transcendental correction away from the immediate neighborhood of the 
turning points f must be linear combinations of 


(83) 


iij- v V" * F *' (' ,,+9 ‘ " '■!"* 

s 


The above asymptotic expressions for two linearly independent solutions to the homoge- 
neous part of (80) simplify considerably, once we realize that the steady state function 
f ~ 7/1 away from the immediate neighborhood of the critical points as pointed out in 
Part 1. Using the expressions for /1 , one finds that to the same order of approximation to 
which (83) is valid, the leading order transcendental correction are multiples of w/f, and 
given by 


(84) 




(ji - p 2 ) 3/s ( i+^v 1 ( 
(i-pV) 3/8 V 1- <*$) U + °7 


where 

(85) 


g(l - aV) 
1 4a(l — a*) 


( 86 ) 

and 

(87) 


, °(» ! - 1?) 

2 4a(l — a*) 


LM _ (1 - pV)!V - ifl 

(<r 2 - a 2 ) z (l - a 2 f 2 ) 2 


as before in the case of the steady state. We notice the similarity of (84) with the expres- 
sions for gi and g-i in part I (Equations 33,34). Each of u> Hl and uji, is transcendentally 
small or large depending on the sign of Re P where, as in part I, 


p(<) = i r d S lw(s) 

Ji/p 


Thus the Stokes lines picture of part I of this paper is critically relevant in determining 
the nature of transcendental correction. Our work in the 1st part of this paper has shown 
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that steady solutions exist only when a 2 — p 7 > 0 . Further, we found that at sufficiently 
small surface tension, U on any branch of solution approaches 2 from below and so in 
the limiting case of small 7 , one can assume that p 2 is positive for fixed a . Thus only 
Fig. 1 in part I of the paper is relevant in the stability analysis in the limit of zero surface 
tension. 

On inclusion of the leading order transcendental correction, the leading order asymp- 
totic behavior of u in sector I (see Fig. 1) is given by 

(89) u ~ Uo + CiUjj l + C%<jj h 3 

As mentioned in section 4 of Part I, Re P has a maximum value of 7 1 / 2 /? in sector I. 
Thus in order that each of the multiples of u;#, and uh 3 be transcendentally small in 
all of sector I, it is necessary to have 

(90) Ct e fi = 0(1) 

(91) C 2 = 0(1) 

In sector II 

(92) u> ~ Wo + Czuh 1 + 

In order that the transcendental corrections in all of sector II be small it is necessary that 

(93) C A e~' 1 ~ 1/3Pm = 0(1) 

(94) 0 3 = 0(1) 

In sector III, the leading order behavior is given by 

(95) oj ~ u>q + 0$u> jj 1 

Note that it is not possible to have a multiple of u>jj l in each of (95) since it grows without 
bounds at f = — a inside the unit f semi-circle. In order that (95) be a transcendentally 
small correction everywhere in sector III of Fig.l, it is necessary that 

(96) C 6 e y/ap "* = 0(1) 

Now, we assume that a — p small with a of order unity, as assumed in Section 4b of 
part I of this paper. As before in part I, section 4b, we now introduce the local variables 
near f = -1/p by defining 

(97) 1 + af = e Xi 
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ea 


(98) /(f) = ; D(x.) 

Then to the leading order in e , equation (80) transforms into 

« M * (i-K) S * is 


where 

( 100 ) 

where 

( 101 ) 


M = (l + xx + Xil)') 172 


/?' = 


7(1 - a 4 ) 3/4 


and the function I>(xi) is determined as described in section 4b of part I. For large 
x\ with corresponding f in sector II of Fig. 1, 

( 102 ) 

w(xi) ~ Xi /(4a) ^ +A 1 (/30 x 7 s/8W{4 “ ) e-^' 1/2 ( 1+s ^ 3/4 / 3 + A 2 ()90^r 3/8W(4a) e 4/3 ' 1/2(1+Il)3/4/s 
and this matches with the expression for u in (92) in sector II provided 

J p A ^/l/2 ^ 

(103) — e"^ = — = e i5ir/4 a~ i/4 ( 1 - a 4 ) 3/ V /(4 " )_3/8 2~ 3/8 2“ a/4a 

C 3 C 4 

where J? is defined as in section 4b, part I (Equation 83). 

For large Xi with corresponding f in sector III, 

(104) 

and this matches with (95) provided 

(105) A 4 = 0 
and 

(106) ^ = e’ 3 ,r / 4 a- 3 / 4 (l - a <) 3 / 8 ^/(4«)-3/8 2 -3/8 2 -a/4a 

C 3 

We now move to the immediate neighborhood of f = - Introduce local variables 


(107) 

(1- 

(108) 

/(f) = 


f (f) = (1 T a 4) D *( x *) 
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where D 2 {x 2 ) is determined as described in section 4b of Part I. To the leading order in 
e , equation (80) becomes 


(109) ^ + (- - ) JT + f 1 ~ | ~PM 2 D 2 + 0M 2 s ) = 0 

dx J \x 2 2 x 2 } dx\ \ 2 ) x\ dx 2 2 a x\ 


For large values x 2 , with corresponding f in sector II, 

(no) 


and this matches to behavior of solution given in (92) in sector II provided 

A.ipfiR A.-ifr/tR 

(111) dil_ = _ a 4 ) 8/8 £ -./(4a)-3/8 2 -3/8 2 ./4 a 

C 3 C 4 

For large values of x 2 with corresponding f in sector I, 

(112) 

and this matches with (89) provided provided 

(113) 21! = = a -»/4(l _ a 4 ) 3/8 £ _ < ,/(4a)-3/8 2 -3/8 2 < T/4« 

Ci c 2 

iFrom (94) and (103), it follows that A\ is transcendentally small compared to A 2 since 
t 1 1/2 Pm is transcendentally small and much smaller that terms involvint /?' since under 
the assumption made here, a — p « 1 . Thus we can neglect A\ altogether. Then 
together with (105), this implies that we are interested in solution to. (99) so that for large 
x\ with corresponding f in sectors II and III in Fig. 1, the transcendental correction 
is small. This would determine a unique solution to (99), as will be shortly argued. As 
far as (109), we note from (93) and (ill) that Aq is transcendentally small compared to 
As and will therefore be set to 0. This means that we should require that the solution to 
(109) contain no transcendentally large term for large x 2 when the corresponding f is 
in sector II. 

We now have to use the condition of smooth bubble back. When a is real, this would 
imply that the mode w(f,o) be real in some open interval containing f = 1 . From 
(80) , it is clear that if w is real on such an open interval, then it must be real for f in 
(a, 1 /a) as all the coefficients of the differential equation are real in that real f interval. 
This implies that the solution u(x 2 , a) of (121) must be real on the postive real x 2 axis. 
However, when a is complex, we cannot demand that each mode u(x 2 ,o) should be 
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real on the real positive x 2 axis. The condition of a smooth bubble back is equivalent to 
requiring that on the positive real x 2 axis 

u(x 2 ,o) = [w(x 2 ,ct*)]* (114) 

On formally taking the complex conjugate of equation (109) on the positive x 2 axis and 
relating the corresponding solution to the condition (114), it is easily argued that the we 
must have the property 

w(x 2 ,a) = [oj{x* 2 , a*)]* (115) 

for any x 2 not necessarily on the positive real axis. Thus the condition that for large 
x 2 corresponding to f in sector II, the solution to (109) does not contain any transcen- 
dentally large corrections implies from (115) that for large x 2 with corresponding f in 
sector VII (which is a reflection of sector II in Fig. l) also contain no trancendental large 
corrections. Indeed, this condition is equivalent to (115) since if we were to require that 
the asymptotic behavior for large x 2 in sectors II and VII contain no transcendentally 
large terms, we will be assured that such solution will automatically satisfy the condition 
(114) and therefore (115) in view of the symmetries of (109). ’ But the condition of no 
transcendentally large terms in two sectors will uniquely determine a solution to (109) as 
argued earlier for the steady problem in part 1. Note that in the special case when a is 
real, such a requirement is indeed equivalent to requiring that u>(x 2 ) be real on the real 
and positive x 2 axis. The actual solution in (109) is not important as it is not needed 
for the leading order determination of the eigenvalue and the corresponding eigenmodes 
as the tip and the back problems decouple to the leading order. 

To ascertain that the tip of the bubble remains smooth under time dependent symmetric 
perturbation, we must require that w(f, r) as defined originally should be real in some 
open interval on the real ? axis containing -1. As before with the neighborhood of 
C = +1 , it is easily argued from (80) that if indeed such an open interval exists then 
w(f, t ) is analytic in the entire interval (—1/a, —a) . However, the same is only true 
for each of the modes w(f , a) if a is real. Generally, for complex a , the condition of 
smooth tip is equivalent to 

u(x u a) = [w(*i, a*)]* (116) 

on the positive real Xi axis, where oj(xi, cr) is the solution to (99) satisfying conditions 
in sectors I and II as mentioned earlier. It is easily argued that the condition (116) is 
equivalent to requiring that there exists solution to (99) with no transcendentally large 
term for large xi with corresponding f in sectors I, II, V and VI , where sectors V and 
VI are the reflections of sectors I and II on the part of the real f axis between a and 
1/a (see Fig. 1). From the figure, it is clear that this means that the solution to (99) 
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satisfying the smooth tip condition is equivalent to requiring that there be contain no 
transcendentally large term for large x x in the entire interval (— 7r, 7r) . Note that such 
a requirement is also valid when a is purely real as can be seen by taking the complex 
conjugate of equation (99). 

We now relate the solution to (99) satisfying the above conditions to a problem which 
we already solved in the context of the linear stability of a finger (Tanveer, 1987c). We 
introduce the transformations 



£ = -is!' 2 

(117) 


D = 2D 

(118) 


II 

(119) 

Then equation (99) becomes 



d 3 u 1. 3 ~ . d 2 u) ( 

3 + £ ^ 2 *^ MD ^ ~d& + ( 

1 - ^ipMD + ip'M 8 ) 

1 dw .a . . 

?dxi + t a £ 3 ^ ~ ° ( 12 °) 


where 

M - 

The requirement of no transcendentally small terms for argument of Xi in the interval 
[ — 7r, 7r] is equivalent to requiring that in the £ plane for large argument, there be no 
transcendentally small terms for argument of £ in the interval [ — tt, 0] . This is precisely 
the problem that was encountered for the finger (Tanveer, 1987c) and numerical solutions 
were found for the cases when /?' was of order unity. Only one branch of solutions for 
which /?' was the smallest, i.e. U the largest was found to be stable and all the other 
branches unstable to tip breaking disturbances. Given the correspondence we have just 
established, the same must be true for the bubble. For large /?' , ( 6 in the notation of 
the Tanveer paper) , asymptotic expressions were found for a and it was found that there 
were many unstable eigenvalues. 

In addition to what has been reported in the Tanveer (1987c) about the limiting values 
of a as 7 — > 0 , we note that in the context of the finger eigen values a for which 
Re a < 0 is not physically acceptable since they correspond to singularities in the tail of 
the finger. In the case of the bubble, these modes physically correspond to disturbances 
affecting the sides and the back of the bubble, as seen before numerically (Tanveer & 
Saffman, 1987b). Using the procedure described before (Tanveer, 1987c), we found two 
nonzero pairs of discrete eigenvalues a = —1.65 ±0.004* and a = — 4.1060±0.11294» for 
the Mclean-Saffman branch. 

4. Discussions and Conclusions 
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Here in this part of the paper, we have demonstrated how the calculation of transcen- 
dentally small terms in surface tension is crucial in the proper prediction of eigen values 
and corresponding eigenmodes for the linear stability of a bubble. We have shown that 
only one branch of bubble solutions is stable, while all the others are unstable, a result 
supported by earlier analytical and numerical work on the finger. The next step will be 
to calculate modal interactions in for calculation of non-linear interactions. It is to be 
pointed out that in this problem, the calculation of stability and the modal interactions in 
an analytical procedure is very important since all the usual numerical procedures break 
down because of inherent ill posedness of the time evolution problem as surface tension 
tends to zero. 
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Figure 1 Stokes Lines in the f plane for a 2 > p 2 > 0, 
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